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The dependencies of the melting point and the lattice parameter of supported metal 
nanoclusters as functions of clusters height are theoretically investigated in the 
framework of the uniform approach. The vacancy mechanism describing the melting 
point and the lattice parameter shifts in nanoclusters with decrease of their size is 
proposed. It is shown that under the high vacuum conditions (p < 10~ 7 torr) the 
essential role in clusters melting point and lattice parameter shifts is played by van 
der Waals forces of cluster-substrate interaction. The proposed model satisfactorily 
accounts for the experimental data. 

The investigation of the properties of supported metal nanoclusters attracts now the increased 
attention, first of all motivated by the need in the fundamental physical basis for the further development of 
nanotechnology related to single-electron and nanoelectronic devices, new electronic and magnetic 
nanomaterials, nanocluster catalysts and nanofilms. 

One of the interesting properties of nanometer-sized systems is a shift of the melting point and the 
lattice parameter with decrease of cluster size (size effect). This effect was experimentally observed for Au, 
Ag and Cu clusters deposited onto various substrates (W, C) [1-6]. The existing theoretical descriptions of 
the observed size effect are based on various approaches [1, 2, 7-11]. However, their common and essential 
drawback is that the influence of the substrate is not taken into account. At the same time it is 
experimentally proved [6] that substrate plays an important role and the use of different substrates (W, C) 
results in the significant difference in melting points of clusters (up to 20% for 30 A clusters on carbon and 
tungsten substrates [6]). Physically it can be related to the additional energy of the supported cluster due to 
the contact potential difference [12], the lattice mismatch at the cluster-substrate interface [13], and the 
cluster-substrate van der Waals interaction [14] that can lead to the shifts of cluster's lattice parameter and 
melting temperature. Therefore, the understanding of the physical reason for the size shift of cluster's 
melting point and lattice parameter needs the determination of the influence of the substrate on cluster's 
properties. 

In this paper the theoretical description of the melting point and the lattice parameter shifts for the 
supported nanoclusters of various metals carried out in the framework of the uniform approach using the 
vacancy mechanism taking into account the cluster-substrate interaction is proposed. In the framework of 
this mechanism the change of cluster properties with decrease of its size is described as a result of 
anharmonic oscillations of cluster atoms, leading to the decrease of the energy of vacancy formation and, 
hence, to the additional creation of vacancies in cluster with decrease of its size. 

The estimations show that the change of the energy of vacancies formation in a cluster with size of 
-30 A due to the van der Waals interaction with substrate can be as much as 8fj.^f^ ter | ft _ 30 ^~0.2eK , 

while the energy change due to the cluster-substrate lattice mismatch, to the contact potential difference at 
the interface, as well as to the adsorption of atoms at cluster surface [16] under the pressure (P < 
10 -7 torr) does not exceed 0.05eV . 

Therefore it can be concluded that the major effect on the shift of the melting point and the lattice 
parameter of supported nanoclusters under the high vacuum conditions (P < 10 -7 torr) is produced by the 
van der Waals cluster-substrate interaction. 

This effect can be described considering the van der Waals interaction of a supported metal cluster 
with bulk metal or semimetal substrate [12,14]. Van der Waals forces are of the electromagnetic character 
and arise due to the mutual polarization of neutral atoms of the interacting objects. The presence of these 
forces can lead to the change of the chemical potential of a cluster. In the case the change is positive, it can 
be interpreted as an additional thermodynamic potential per atom of a cluster, arising due to the interaction 
with substrate. It results in the decrease of the binding energy of atoms in a lattice and, therefore, in the 
increase of the probability of vacancies formation in a cluster. Besides van der Waals forces, the formation 
of vacancies in a cluster occurs due to the presence of cluster surface, that results in the dependence of the 
energy of vacancy formation on its distance from cluster surface, as well as on cluster shape [17]. Under 
assumption that the clusters under investigation [1-5] satisfy the requirement l/h > 1 (where 1=0. 7-10.0 
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nm and h=0. 3-5.0 nm are cluster's lateral size and height, consequently), the cluster surface can be 
considered as a flat and the influence of the boundary effects related to cluster's shape can be neglected. 

The equilibrium concentration of vacancies in a cluster can be found from the requirement of the 
minimum of the free energy of cluster's vacancies subsystem taking into account the interaction with 
substrate: 

r =° <" 

on n=n 

where n is the average concentration of vacancies in a cluster. The interaction of vacancies at the distance 
more than the character size of the vacancy (the possibility of divacancy formation) is attractive, while at 
the distance less than the characteristic size of the vacancy it has a character of a solid core (it is impossible 
to locate two vacancies at the same site of the lattice). Assuming the equilibrium concentration of vacancies 
much less than the atomic density in a cluster, the vacancies subsystem of a cluster can be considered as a 
van der Waals' gas [18]. In this case the free energy of the vacancies subsystem after decomposition for 
virial coefficients is given by [19]: 

F = F W +^;S(T) = i/(l-exp(-^))^ (2) 

where TV is the number of vacancies in a cluster, V is the cluster volume, U l2 is the interaction potential of 
vacancies, and F id is the free energy of non-interacting vacancies system. Considering U 12 as a rectangular 
potential well of depth E d (the energy of divacancy formation), the virial coefficient B (T) is expressed as: 

B(T) = b + 7b(l-exp(^-)) (3) 

where b = ^p- , r is the characteristic size of a vacancy. The free energy F id of non-interacting vacancies 
system is given by [17]: 

Fid = / E v ndV + Tf [n In Q - n] dV ; E v = E* - 8pL ; (4) 

E * = h J ( Eb ~ (<Eb ~ Es ^ 6XP (~ 2a)) dX 
Where E B ,E S are the energies of vacancy formation in the volume and at the surface of a cluster, 
consequently, n is the atomic density for a cluster, a is the cluster lattice constant, Sfj. is the additive 
chemical potential of a cluster arising due to its interaction with substrate [14]. In the expression (4) the 
energy of the vacancy formation is assumed depending on its distance from cluster surface, and E„ is the 
energy of vacancy formation averaged for a cluster. Substituting the expression (4) into (2), the free energy 
is obtained: 

F = f E v ndV + Tf n In (^) - n] dV - nT ln(l - nb) - 7n 2 bT (l - exp (^ffj (5) 
The equation for the equilibrium concentration of vacancies n in a cluster follows from (5) and (1) as: 

E* - 8 ii + T\n^- - Tln(l - nb) + nT-^- lAnbT (l - exp (^ffj = (6) 

Where Sfi can be calculated if the complex dielectric constants of substrate and cluster are known [14,12]. 
The effect of the cluster surface taken into account in the expression for the energy of vacancy formation, 
the value <5/i can be calculated for a film of thickness h deposited onto substrate [14,12]: 

_ a 3 7uo __ .re (gfOQ-l^gfOQ-gsOO) 

d v ~ Z 8***3 ; " - Jo ( £/ao+l) ( £/ao+£sao ) d $ W 

Where if = a) is the imaginary part of the electromagnetic field, s s , £t are the dielectric constants of 
substrate and cluster materials, consequently, and z is the number of the nearest neighbors in a cluster. 
Under assumption of the validity of the expression Ej (if) = 1 + — p-,y = s,f for the dielectric constants 
Sf, s s [12,20], the value <5/i is given by: 

r z Ann 6 o s \ Of ) y 
where Of, a s are the conductivity of cluster and substrate, consequently. 
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The conductivity o s of a cluster can be calculated by taking into account the contributions to the 
intrinsic resistivity of the cluster material p and the additional resistivity p vac arising due to the additional 
scattering of conduction electrons from vacancies. The additional resistivity p vac depends on the 
concentration of vacancies defined by the additional chemical potential arising due to the presence of a 
substrate. Thus, the expression for the cluster conductivity takes the form: 

Of = (Po + Pvac)' 1 I Pvac = Pvac ^ (9) 

Where p^, ac is the additional resistivity per a vacancy that does not depend on n. The expressions (6), (8) 
and (9) allow to calculate the equilibrium concentration of vacancies in a cluster n taking into account the 
cluster-substrate interaction. The theoretical dependence of the relative concentration of vacancies for Au 
cluster on graphite surface as function of cluster height at room temperature (7=300 K) calculated from the 
equations (6), (8) and (9) is plotted in Fig. 1. 
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Fig 1. The equilibrium concentration of vacancies n in Au cluster on graphite at 7=300 K relative to the atomic density n 

calculated as function of cluster height h. 



It is seen that the concentration of vacancies sharply increases for cluster height less than 22 A, that 
can result in the modification of the modulus of elasticity of cluster and, as a consequence, to the shifts of 
the lattice parameter and the melting point [18]. 

The melting point can be calculated taking into account the dependence of the modulus of elasticity 
on the concentration of vacancies in a cluster [18]. Using the vacancy model of melting [14] and 
considering vacancies as dilatation centers, it is possible to calculate the shear modulus K' as function of the 
concentration of vacancies [14]: 



~fl-—l 

K' =K-^n\KAV) 2 " l " n) 



(10) 



Here K is the shear modulus of the defect free material and AV = 0.1 h- 0.2a 3 is the dilatation volume. 
Defining the melting point of a cluster from the Born's criterion K (T melt ) = [18], the dependence of the 
melting point on cluster height can be obtained from the equation (10). The dependence of the melting point 
of Au clusters supported on carbon and tungsten substrates as function of cluster height, calculated from (6), 
(8), (9), (10), as well as the dependence calculated without the influence of the substrate, are presented in 
Fig. 2. 
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Fig. 2. The dependence of the melting point T meit of Au cluster on its heigh h. Solid and dashed curves are the theoretical 
dependences calculated using equations (6), (8), (9) and (10) for Au clusters on W and C, consequently. Dotted line is the 
dependence calculated without the influence of a substrate. Points are the experimental data [1,6]. 

The plots in Fig.2 evidence that the substrate plays an important role in the shift of the melting point 
of clusters. It is seen that different substrates result in the significant shift of the melting point of clusters 
that is satisfactorily described by the proposed model (see Fig.2). 

The formation of extra number of vacancies in a cluster compared to the bulk metal leads to the 
reduction of the lattice parameter of a cluster that is related to the vacancy concentration [21]: 



— = - — n (11) 

a 3Vb v ' 

Where a is the lattice parameter in the absence of vacancies, Aa is the change of the lattice parameter due 
to vacancies, V is the atomic volume of a defect free crystal, and V x is the change of the atomic volume 
caused by the presence of one vacancy. Using the relationship between V and n the lattice parameter in 
the presence of vacancies can be expressed as: 

(12) 



a _ ^ In 

a 3 n 



The expression (12) shows that the additional formation of vacancies in a cluster results in the decrease of 
the lattice parameter. 

Equations (6), (8), (9), (10) and (12) allow to calculate the dependence of the lattice parameter on 
cluster height. Such dependences for the clusters of Cu, Au and Ag deposited on graphite at room 
temperature are presented in Fig. 3. The parameters used in the calculations of the dependencies plotted in 
Figs. 1-3 are summarized in Table 1. 
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Fig. 3. The lattice parameter a for Au, Ag (a) and Cu (b) clusters on graphite as function of cluster height h at room temperature. 
Solid lines are the theoretical dependences calculated using equations (6), (8), (9), (10) and (12). Points are the experimental data 

[2]. 



Table 1. Parameters used in the calculations. 



Parameter 


Au 


Ag 


Cu 


The energy of vacancy formation in the volume, eV 


0.9 


1.11 [I5J 


1.28 L1SJ 


The energy of vacancy formation at the surface, eV 


0.25 [23J 


0.16 U5J 


0.35 L1SJ 


Divacancy formation energy [23J , eV 


0.1 


0.1 


0.05 


Lattice parameter [24J , A 


4.08 


4.09 


3.61 


Specific conductivity [23J , 10" 6 Q-cm 


2 


1.5 


1.55 



Fig. 3 shows that the additional formation of vacancies in a cluster due to its interaction with 
substrate results in the decrease of the lattice parameter that agrees with the experimental data. 

The additional energy of a cluster per one atom due to its interaction with substrate can be estimated 
in the frames of several mechanisms taking into account the contact potential difference, the lattice 
mismatch at the interface, the van der Waals interaction with the substrate, as well as the adsorption and the 
chemical reactions at cluster surface. The estimations will be carried out for the cluster of 30 A in height. 
According to the equations (8) and (9), the additional energy of a cluster per one atom due to the van der 
Waals interaction with substrate is 5^f^ ter \ h _ 30 ^~0.2eV. 

The shift of the cluster chemical potential due to the presence of the cluster-substrate interface is 

given by [12]: Sfj. = A exp(— kx) ; A = K2< ^ 2 ; k\ 2 = 47re 2 ^ where \i\,\i2 are the chemical 

potentials of graphite (1) and Au (2), n° 2 is the electron density, k = — is the reverse Debye screening 
radius. 

The additional energy of a cluster obtained due to the contact potential difference can be estimated 
assuming that the influence of a substrate is considerable only within the first monolayer of atoms at the 

interface: Sfj. c ~^ f SfidV ~ 5 ^ 2 J K > where V, h are the volume and the height of a cluster. In this 
case the value of the additional energy due to the contact potential difference equals to 

Cluster UoA - -°- 05 ^ 
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The additional energy of a cluster per one atom due to the lattice mismatch at the interface can be 
estimated considering that the characteristic temperature of the transition to the disproportionate phase is 
about 500 K. [26]. Thus, the additional energy per one atom due to the this mechanism is 

The maximal additional energy per one atom due to the adsorption and the chemical reactions at the 
cluster surface can be estimated as <5/i^ ter = N t Ef r , where N t is the number of atoms adsorbed at 
cluster surface, and Ef r is the hemisorption cohesive energy (Ef r ~l h- l.SeV). Fractional coverage of the 
cluster surface is 0.001 under P = 10~ 9 torr , so the additional energy per one atom of cluster due to the 
adsorption and the chemical reactions is less than 0.01 eV. Note that this value increases with pressure and 
under P = 10~ 5 torrit may reach, ~0.7eV that may be the reason of the experimentally observed [2,27] 
lattice extension under the low vacuum conditions. 

Summarizing the above estimations, it is seen that under the high vacuum conditions 

Spinster | /=3oA » 5 ^dus^r • s Ater • s Vciuster • Therefore under these conditions in the framework of the 
considered mechanism the shift of the melting point and the lattice parameter of clusters with decrease of 
their size is defined by the van der Waals cluster-substrate interaction. 

The work has been partially supported by the Russian Foundation for Basic Research, grants No. 
02-00759a and 07-02-0 1372a. 
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